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The properties of butane (C4H10) lubricants confined between two approaching solids are inves- 
tigated by a model that accounts for the curvature and elastic properties of the solid surfaces. We 
consider the linear n-butane and the branched iso-butane. For the linear molecule, well defined 
molecular layers develop in the lubricant film when the width is of the order of a few atomic diam- 
eters. The branched iso-butane forms more disordered structures which permit it to stay liquid-like 
at smaller surface separations. During squeezing the solvation forces show oscillations corresponding 
to the width of a molecule. At low speeds (< 0.1 m/s) the last layers of iso-butane are squeezed 
out before those of n-butane. Since the (interfacial) squeezing velocity in most practical applica- 
tions is very low when the lubricant layer has molecular thickness, one expects n-butane to be a 
better boundary lubricant than iso-butane. N-butane possessing a slightly lower viscosity at high 
pressures, our result refutes the view that squeeze out should be harder for higher viscosities, on 
the other hand our results are consistent with wear experiments in which n-butane were shown to 
protect steel surfaces better than iso-butane. 



I. INTRODUCTION 

Modern materials are subject to increasing loads and 
are used under still more demanding conditions. This 
emphasises the need for a better understanding of fric- 
tion, lubrication and wear phenomena 0,0. 

An example of the above is the fuel lubricated diesel 
engine injection pump, capable of delivering in excess of 
2000 bars today compared to 500 bars 15 years ago. This 
increase in injection pressure has ensured higher engine 
efficiencies and lower pollution levels. 

Historically the diesel oil sulphur reduction both in the 
1960s and 1990s combined with the high pressures have 
played a significant role in injection pump durability. Sul- 
phur, along with other polarity inducing atoms, is part of 
polar species that ensure low wear in boundary lubrica- 
tion, a predominant regime in pumps lubricated by diesel 
oil. Today there exist several accelerated laboratory tests 
capable of predicting the lubricating abilities of a diesel 
oil. If a fuel fails such a test, lubricity additives in small 
proportions are added to ensure proper wear resistance. 

The currently used diesel oil wear test in Europe is the 
high frequency reciprocating rig (HFRR) which is cov- 
ered by several standards |3, |^]. The principle in the 
HFRR is the sliding of a fixed steel ball, loaded with 2N 
on a steel disk for 75 minutes. The motion is reciprocat- 
ing with a stroke of 1 mm and a frequency of 50 Hz. As 
the specimen contact is fuel lubricated the resulting wear 
scar diameter on the ball expresses the lubricity of the 
tested fuel. The configuration of the HFRR is shown in 
figure n 

The HFRR ability of predicting the lifetime of injec- 
tion pumps has been questioned when low- viscosity fuels 
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FIG. 1: The ball on disk configuration in a HFRR wear test. 



are tested. It appears that a decrease in fuel viscosity 
requires an increase in lubricity (smaller wear scar) to 
ensure full lifetime of the pumps. This relation between 
viscosity and lubricity was confirmed as a new fuel ap- 
peared in the 1990s: Dimethyl ether (DME). This fuel 
has a viscosity 20 times lower than that of diesel oil and 
it was established that an adequate lubricity level accord- 
ing to the HFRR standards was not sufficient to protect 
the pump surfaces ^, ^ . 

A study 8] has argued that that fluid viscosity is only 
a secondary property in the HFRR. Based on Molecular 
Dynamics calculations it was discovered that the squeeze- 
out of linear alkanes from surface contacts is primarily a 
function of the length of the molecule. Experimental data 
shown in figure [3 demonstrates that branched alkanes 
perform worse in the HFRR than their linear isomers 
although their viscosities are almost the same @. This 
observation confirms that the bulk viscosity may be a 
secondary property in boundary lubrication. 

The friction and wear properties in boundary lubri- 
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FIG. 2: The wear scar diameter for various alkanes as a func- 
tion of the molecular length. The squares are from 1^ and 
the circles from 8] . a) , c) are for linear alkanes and b) , d) for 
branched isomers. 




FIG. 3: The simulated sample: two elastic walls with lubri- 
cant in between are pushed together. The sinusoidal profile 
of the upper surface generates a line contact in the middle, 
where lubricant forms a neck and eventually it is squeezed out 
by the load. Periodic boundary conditions are applied along 
X and y. 
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cation of lubricants with branched and linear molecules 
have been investigated both experimentally 9, 10, 11, 12] 
and usins; theoretical simulation methods [il El El 111 

Whether linear alkanes or their branched isomers 
are the best lubricants according to the cited litera- 
ture remains uncertain. A number of studies indicate 
that n- alkanes are the best lubricants 0, Q, Q 
whereas others claim that branched alkanes perform bet- 
ter 0,0,0,^^. One major weakness of the published 
studies is that the lubricant squeezing is stopped at a sep- 
aration of about eight Angstrom between the surfaces. At 
this point there are still a few monolayers left protecting 
the surfaces, and significant wear originating from cold 
welding should onl y appear when the last lubricant layer 
is squeezed out p^ . 

In the present computer simulation study we investi- 
gate theoretically the lubrication abilities of n- and iso- 
butane extending the squeezing down to a surface sepa- 
ration of zero Angstrom so that the important expulsion 
of the last monolayer is addressed. 

II. THE MODEL 

Figure [31 shows the geometry for the simulation of the 
squeeze-out process. Lubricant molecules are confined 
between two elastic walls. The lower wall, or substrate, 
is flat, while the upper one, called block, has a sinusoidal 
profile, so that a contact region appears in the centre of 
the sample. The liquid lubricant adheres to the walls, 
forming a neck in the central region. When the block 
is pushed towards the substrate, the lubricant has space 
to expand laterally into the vapour. Periodic boundary 
conditions are applied along the x and y directions. The 



FIG. 4: The confined regime: the liquid lubricant get ordered 
into 2D layers when it is confined between atomically flat 
surfaces. The picture shows the side view of iso-butane in the 
contact region between the two walls. Light balls represent 
CH and CH3 groups. 

periodically repeated cell forms a rectangle Lx ^ Ly, with 
Lx = 416 A and Ly = 52 A. The contact region is not 
circular, but forms an infinite long strip parallel to the y 
axis. 

The lubricant molecules are described through the Op- 
timized Potential for Liquid Simulation (OPLS) 20, 21]; 
this potential is known to provide density and viscosity 
of hydrocarbons close to the experimental one. 

Each butane molecule comprises four units (particles), 
each particle corresponding to one chemical group CH3, 
CH2 or CH {United Atom model). The interaction be- 
tween particles of different molecules is described by 
Lennard- Jones potentials. The intramolecular interac- 
tions include two body forces that keep the bond length 
C-C close to 1.53 Angstrom, three body forces impos- 
ing a preferred angle of 115 degrees between the carbon 
atoms, and four body forces favouring a well defined tor- 
sion of the molecules. The four body forces apply to 
the sequence of carbon atoms C-C-C-C, thus they are 
present only for n-butane. 

For the iso-butane molecules we introduced also an anti 
umbrella inversion potential, analogously to what already 
done by Mondello and Grest ^2^ for branched hydrocar- 
bons; the only difference here being that in the iso-butane 
molecule it is not possible to distinguish which carbon 
atom constitutes the lateral branch and which constitute 
the backbone of the molecule. Therefore we symmetrized 
the anti inversion potential. 

The interaction between the beads of the lubricant 



molecule and the walls' atoms is also given by Lennard- 
Jones forces, whose energy parameter eo = 18.6 meV, as 
in ref. [2^, yielding an adsorption energy of about 0.3 eV 
per molecule, comparable to that of butane on Au(lll) 



The elastic energy due to the walls' deformation is also 
involved in the nucleation process that triggers the tran- 
sition from n to n — 1 monolayers of lubricant in the 
confined regime . In principle we can achieve a realis- 
tic description of the walls' elastic response by simulating 
thick walls comprising many layers. Since we necessitate 
proper elastic behaviour up to wavelengths comparable 
with the size of the contact area the walls' thickness must 
not be smaller that the lateral size of the box, Ly = 52 A. 
This indeed would imply too many atoms, slowing down 
the calculations. Instead we adopted the model described 
in ref. |2^: for each wall only the outermost layer of 
atoms is considered; these atoms are connected to a rigid 
surface with springs that take into account both the com- 
pressibility and the shear rigidity of the wall. Similar 
springs connect together the square grid of atoms of the 
wall. The advantage of this approach is double: with 
a relatively small number of atoms it is possible to de- 
scribe the long range elasticity of the walls, and to impose 
a curved profile simply by using a curved rigid surface. 

The rigid surface connected to the substrate's atoms is 
flat and its position is fixed. The upper rigid surface has 
the profile 
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where the amplitude of corrugation A is 20, 40, or 60 
A. The case A = 20 A guarantees a larger contact re- 
gion, but it does not allow enough empty space (or rather 
vapour space region) for the squeeze-out of the last two 
layers; it was used only for comparison with the other 
cases, to ensure that the size of the contact region is large 
enough for the nucleation of the squeeze-out process. 

The springs connecting the walls' atoms to the rigid 
surfaces simulate the elastic response of a gold film of 
thickness 50 Angstrom, comparable with the lateral size 
Ly of the sample. 

The substrate consists of 144 x 18 atoms in a square 
lattice with lattice spacing 2.889 Angstrom; the block 
layer is made of 160 x 20 atoms, with lattice spacing 2.6 
Angstrom. The different lattice spacing is to reduce the 
commensurability between the walls. Both substrate's 
and block's atoms have the same mass of gold: 197 a.m.u. 
In most of the simulations we employed 2800 lubricant 
molecules, which deposit on the two walls forming about 
two monolayers of adsorbate on each surface. 
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FIG. 5: Squeeze-out of butane at T = 300 K. The pressure 
drops when the number of layers switches from 3 to 2 and 
from 2 to 1. Squeezing speed 1 m/s. 
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FIG. 6: Top view of the lubricant molecules in the central 
part of the contact area immediately before the expulsion 
of the last monolayer (top), after 0.5 nanoseconds (middle), 
and after 1 nanosecond. T — 300 K, squeezing speed 1 m/s. 
(Block and substrate are parallel to the plane of the figure) 



III. SIMULATION RESULTS 

Two simulations are prepared in parallel for iso- and n- 
butane with identical conditions. Initially the whole sys- 
tem is thermalized and the lubricant adheres to the two 



walls. Then the rigid upper surface is moved down to- 
wards the substrate with constant speed, i.e., zq in equa- 
tion is constrained to decrease linearly with time. As 
the lubricant layers on the two walls get in contact, a 
neck is formed in the contact area. The lubricant clearly 
shows layering both for the n-butane and for the more 



time (ns) 

FIG. 7: Pressure versus time during the complete squeeze- 
out of three monolayers at T = 350 K, squeezing speed 1 m/s. 
The pictures show the lubricant while the last monolayer is 
being removed. Corrugation of the upper profile: A = 60 A; 
contact area: 10% of the cell's area. 



disordered iso-butane (Fig. The change in the num- 
ber of monolayers takes place abruptly, with a sudden 
decrease of pressure due to the relaxation of the elastic 
energy stored in the walls. 

Fig.[51shows the behavior of pressure versus time while 
squeezing at speed 1 m/s and at temperature 300 Kelvin. 
The y-axis contains a spatial average of the pressure, that 
is the vertical force divided by the area LxLy of the sim- 
ulated cell. The contact region is about 15% of the whole 
cell size, and the real pressure in the center of the con- 
tact area is much larger than the average pressure. Nega- 
tive pressures observed in some time intervals are due to 
the attractive capillarity forces of the neck of lubricant. 
The plotted curves show that iso-butane is squeezed out 
earlier and at a lower pressure. This is the typical be- 
haviour we observed when there is more than one mono- 
layer. Conversely the removal of the very last monolayer 
appears to be easier for n-butane. As can be observed in 
Fig. ini the ejection of the iso-butane not only happens 
later, but it is slower and a small group of 7 molecules 
remains trapped at the end. This suggests that the sim- 
ulated squeeze-out of the last layer is strongly affected 
by dynamical effects, such as the drag friction when the 
lubricant slides on the solid walls. Thus, even though 
the nucleation barrier of the hole may be smaller for 
iso-butane, the simulated system might not have enough 
time to exit its metastable state before the pressure is 
further increased by the vertical (squeezing) motion of 
the block. 

The behaviour at higher temperatures is qualitatively 
the same. At T = 350 Kelvin the switch from 3 to 2 






FIG. 8: Side and top views of the iso-butane monolayer at 
different time (and pressure). Images refer to the simulation 
of Figure 



monolayers and from 2 to 1 monolayer happens first for 
iso-butane, but the very last layer of iso-butane stays 
longer between the walls, and 5 iso-butane molecules re- 
main trapped there (see Fig. [TJ. Another independent 
simulation at T = 330 K (not shown) confirms this trend. 

The transition from 2 to 1 iso-butane monolayer reveals 
another interesting feature: the pressure drops down in 
two stages, as is shown in the plot of figure [3 for time 
close to 1.2 ns. What happens is clearly illustrated in 
the snapshots of figure [SI initially there is the squeeze- 
out of one of the two layers. The remaining monolayer 
does not have a preferential orientation of the molecules 
parallel to the walls, the plateau in the pressure versus 
time graph at t ~ 1.2 ns is due to this thick monolayer. 
Finally the molecules change their orientation forming 
a thinner monolayer, yielding a smaller second pressure 
drop. This kind of two-stage squeezing of the second layer 
has been observed in other simulations too, but only for 
iso-butane. It is in fact a steric effect due to the shape 
of the molecule. Similar pressure-induced phase transi- 
tions have been observed in other computer simulations 
of squeeze-out. Thus, in Ref. 27] it was shown that Xe 
atoms between two solid surfaces exhibited triangular bi- 
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FIG. 9: Last monolayer squeeze-out at low speed. T = 300 
K, squeezing speed 0.03 m/s. Iso-butane is removed earlier at 
a lower pressure, but the cluster of six molecules shown in the 
picture remains trapped up to t = 3.6 ns. Below: the average 
position of the upper wall's atoms versus time. 
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FIG. 10: Pressure required to squeeze-out the last lubricant 
layer while pushing at speed 0.03 m/s. Data refers to the 
spatial average over the whole simulation cell. The pressure 
in the centre of the contact area is much larger. Each point is 
the result of a single independent simulation, requiring about 
2 mounths of CPU-time on a standard PC. 



layers at low squeezing pressure, which abruptly trans- 
formed to fcc(lOO) layers parallel to the solid surfaces 
when the pressure increased. This transition allowed the 
surfaces to move closer to each other which released elas- 
tic energy. 

The simulations at lower squeezing speed give clear 
evidence that the the relatively higher lubricity of the last 
layer of iso-butane is a kinetic effect. Thus, at T = 300 
K and for the compression speed 1 m/s (Fig. 13 n-butane 



is removed first, while when the speed is 0.1 m/s the 
removal of the two lubricants happens almost at the same 
time (not shown). Finally when the squeezing is carried 
at the speed 0.03 m/s iso-butane is removed first and 
at a lower pressure (Fig. In this case there is still 
some small island of molecules trapped, shown in the 
picture, confirming once more that the lateral sliding of 
iso-butane is slower. 

A further evidence that the removal of the last mono- 
layer is influenced by kinetic effects comes from the 
squeezing pressure versus pushing speed Vz'. keeping all 
the other conditions identical, we observed that the pres- 
sure at which the last layer is removed does decrease when 
the block is pushed down slowly. For example, for iso- 
butane at T = 300 K we obtained the following pressures: 
0.340 GPa when = I m/s, 0.257 GPa when = 0.1 
m/s, and 0.239 GPa when Vz = 0.03 m/s (as usual, these 
are spatial averages, the true pressures in the contact area 
being much higher). Similarly, the corresponding simu- 
lations for n-butane gave 0.321, 0.252, and 0.263 GPa 
respectively. The lubricant is able to nucleate the hole 
for the squeeze-out at the lowest of these pressures, but 
at high pushing speed there is not enough time to exploit 
this possibility before the pressure is further increased. 

It is interesting to observe what happens when the tem- 
perature is changed, particularly considering that in re- 
ality the microcontacts where the squeeze-out happens 
are likely to be much warmer than the environment, at 
least during sliding. The temperature is likely to speed 
up the mobility of the molecules, reducing the influence 
of the drag force of the walls. Moreover the larger ther- 
mal fluctuations should strongly favour the nucleation 
of the critical hole in the lubricant layer. Both effects 
tend to reduce the pressure needed to squeeze the lubri- 
cant. This is indeed the trend which we found and that 
is observed in Fig. On the other hand, any difference 
AG between n-butane and iso-butane in the free energy 
barriers for nucleation of squeeze-out should become less 
important as the temperature increases, as the ratio be- 
tween the nucleation probabilities is mainly affected by 
the Boltzmann factor exp(— /3AG), which goes to 1 when 
T ^ oo. 

Unfortunately the low temperature conditions are hard 
to analyse through simulations: the randomness of the 
single squeeze-out event is enhanced by the strong de- 
crease of the Boltzmann factor. Actually the transition 
time for a nucleating process has an exponential distribu- 
tion, where the standard deviation is equal to the average 
value. Thereafter the data of a single simulations become 
no more reliable and only an average over many indepen- 
dent simulations — not easily affordable due to computer 
time limitations — would show the real trend. We do not 
have enough statistics in Fig. for a reliable compari- 
son of the behaviours of n-butane and iso-butane below 
room temperature. 

Finally we consider the effect of lateral sliding of the 
walls with respect to each other. We ran some simula- 
tions (not shown) at different temperatures with squeez- 
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ing speed and lateral speed both equal to 1 m/s, and 
we observed that the removal of the last monolayer hap- 
pens almost simultaneously at the same pressure. The 
relative motion of the walls tends to take away the ef- 
fect of the drag between walls and lubricant, which is 
indeed responsible of the higher lubricity of iso-butane 
when squeezing fast. This is an important remark, since 
in many practical situations there is lateral sliding. 

The significance of the temperature for the squeeze- 
out of the two different butanes is of outermost impor- 
tance for future wear experiments. The results in figure 
121 are all obtained at 298K. At this temperature and at 
the given sliding conditions it seems that iso-butane is 
more easily squeezed out from the ball-disk contact than 
n-butane as the resulting wear is slightly lower for the lat- 
ter one. In real injection pumps the temperature could 
be very low (at start up) and very high (running condi- 
tions). The present simulation results imply that wear 
protecting properties of different molecules in the fuel 
may not react the same way to changes in temperature 
and sliding conditions. This means that a wear test may 
not reflect reality unless test conditions such as temper- 
ature and sliding velocity are varied. This new aspect 
of wear experiments should be looked into in the near 
future. 



IV. SUMMARY AND CONCLUSIONS 

We have studied the properties of butane (C4H10) lu- 
bricants confined between two approaching solids using a 
model that accounts for the curvature and elastic prop- 
erties of the solid surfaces. We considered both linear 
n-butane and the branched iso-butane. In the case of the 
linear molecule well defined molecular layers develop in 
the lubricant film when the width is of the order of a few 



atomic diameters. The branched iso-butane forms more 
disordered structures which permit it to stay liquid-like 
at smaller surface separations. 

During squeezing the solvation forces show oscillations 
corresponding to the width of a molecule. At low speeds 
(< 0.1 m/s) the last layers of iso-butane are squeezed out 
before those of n-butane. Since the (interfacial) squeez- 
ing velocity in most practical applications is very low 
when the lubricant layer has molecular thickness, one 
expect n-butane to be a better boundary lubricant than 
iso-butane. This is consistent with wear experiments in 
which n-butane was shown to protect steel surfaces better 
than iso-butane. On the other hand, n-butane possess- 
ing lower viscosity at high pressures, our result refutes 
the view that squeeze out should be harder for higher 
viscosities. At high squeezing velocity, the squeeze-out 
of the last monolayer of iso-butane occurs at higher pres- 
sures than for n-butane. We interpret this as a kinetic 
effect resulting from the lateral corrugation barrier expe- 
rienced by the molecules. As the temperature increases, 
the squeeze-out occurs at lower applied pressure for both 
n-butane and iso-butane. This is the expected result 
based on the nucleation theory of squeeze-out 
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